Molecular mapping of the broad bean wilt virus 2 resistance locus bwvr in Capsicum annuum using BSR-seq.

KEY MESSAGE: Bulked segregant RNA seq of pools of pepper accessions that are susceptible or resistant to Broad bean wilt virus 2 identifies a gene that might confer resistance to this devastating pathogen. The single-stranded positive-sense RNA virus Broad bean wilt virus 2 (BBWV2) causes substantial damage to pepper (Capsicum annuum) cultivation. Here, we describe mapping the BBWV2 resistance locus bwvr using a F7:8 recombinant inbred line (RIL) population constructed by crossing the BBWV2-resistant pepper accession 'SNU-C' with the susceptible pepper accession 'ECW30R.' All F1 plants infected with the BBWV2 strain PAP1 were susceptible to the virus, and the RIL population showed a 1:1 ratio of resistance to susceptibility, indicating that this trait is controlled by a single recessive gene. To map bwvr, we performed bulked segregant RNA-seq (BSR-seq). We sequenced pools of resistant and susceptible lines from the RILs and aligned the reads to the high-quality 'Dempsey' reference genome to identify variants between the pools. This analysis identified 519,887 variants and selected the region from 245.9-250.8 Mb of the Dempsey reference genome as the quantitative trait locus region for bwvr. To finely map bwvr, we used newly designed high-resolution melting (HRM) and Kompetitive allele specific PCR (KASP) markers based on variants obtained from the BSR-seq reads and the PepperSNP16K array. Comparative analysis identified 11 SNU-C-specific SNPs within the bwvr locus. Using markers derived from these variants, we mapped the candidate bwvr locus to the region from 246.833-246.949 kb. SNU-C-specific variants clustered near DEM.v1.00035533 within the bwvr locus. DEM.v1.00035533 encodes the nitrate transporter NPF1.2 and contains a SNP within its 5' untranslated region. The bwvr locus, which contains four genes including DEM.v1.00035533, could represent a valuable resource for global pepper breeding programs.

The clinical value of optical genome mapping in the rapid characterization of RB1 duplication and 15q23q24.2 triplication, for more appropriate prenatal genetic counselling.

BACKGROUND: Despite recent advances in prenatal genetic diagnosis, medical geneticists still face considerable difficulty in interpreting the clinical outcome of copy-number-variant duplications and defining the mechanisms underlying the formation of certain chromosomal rearrangements. Optical genome mapping (OGM) is an emerging cytogenomic tool with proved ability to identify the full spectrum of cytogenetic aberrations. METHODS: Here, we report on the use of OGM in a prenatal diagnosis setting. Detailed breakpoint mapping was used to determine the relative orientations of triplicated and duplicated segments in two unrelated foetuses harbouring chromosomal aberrations: a de novo 15q23q24.2 triplication and a paternally inherited 13q14.2 duplication that overlapped partially with the RB1 gene. RESULTS: OGM enabled us to suggest a plausible mechanism for the triplication and confirmed that the RB1 duplication was direct oriented and in tandem. This enabled us to predict the pathogenic consequences, refine the prognosis and adapt the follow-up and familial screening appropriately. CONCLUSION: Along with an increase in diagnostic rates, OGM can rapidly highlight genotype-phenotype correlations, improve genetic counselling and significantly influence prenatal management.

Identification of superior and rare haplotypes to optimize branch number in soybean.

KEY MESSAGE: Using the integrated approach in the present study, we identified eleven significant SNPs, seven stable QTLs and 20 candidate genes associated with branch number in soybean. Branch number is a key yield-related quantitative trait that directly affects the number of pods and seeds per soybean plant. In this study, an integrated approach with a genome-wide association study (GWAS) and haplotype and candidate gene analyses was used to determine the detailed genetic basis of branch number across a diverse set of soybean accessions. The GWAS revealed a total of eleven SNPs significantly associated with branch number across three environments using the five GWAS models. Based on the consistency of the SNP detection in multiple GWAS models and environments, seven genomic regions within the physical distance of ± 202.4 kb were delineated as stable QTLs. Of these QTLs, six QTLs were novel, viz., qBN7, qBN13, qBN16, qBN18, qBN19 and qBN20, whereas the remaining one, viz., qBN12, has been previously reported. Moreover, 11 haplotype blocks, viz., Hap4, Hap7, Hap12, Hap13A, Hap13B, Hap16, Hap17, Hap18, Hap19A, Hap19B and Hap20, were identified on nine different chromosomes. Haplotype allele number across the identified haplotype blocks varies from two to five, and different branch number phenotype is regulated by these alleles ranging from the lowest to highest through intermediate branching. Furthermore, 20 genes were identified underlying the genomic region of ± 202.4 kb of the identified SNPs as putative candidates; and six of them showed significant differential expression patterns among the soybean cultivars possessing contrasting branch number, which might be the potential candidates regulating branch number in soybean. The findings of this study can assist the soybean breeding programs for developing cultivars with desirable branch numbers.

Integrated meta-analysis and transcriptomics pinpoint genomic loci and novel candidate genes associated with submergence tolerance in rice.

BACKGROUND: Due to rising costs, water shortages, and labour shortages, farmers across the globe now prefer a direct seeding approach. However, submergence stress remains a major bottleneck limiting the success of this approach in rice cultivation. The merger of accumulated rice genetic resources provides an opportunity to detect key genomic loci and candidate genes that influence the flooding tolerance of rice. RESULTS: In the present study, a whole-genome meta-analysis was conducted on 120 quantitative trait loci (QTL) obtained from 16 independent QTL studies reported from 2004 to 2023. These QTL were confined to 18 meta-QTL (MQTL), and ten MQTL were successfully validated by independent genome-wide association studies from diverse natural populations. The mean confidence interval (CI) of the identified MQTL was 3.44 times narrower than the mean CI of the initial QTL. Moreover, four core MQTL loci with genetic distance less than 2 cM were obtained. By combining differentially expressed genes (DEG) from two transcriptome datasets with 858 candidate genes identified in the core MQTL regions, we found 38 common differentially expressed candidate genes (DECGs). In silico expression analysis of these DECGs led to the identification of 21 genes with high expression in embryo and coleoptile under submerged conditions. These DECGs encode proteins with known functions involved in submergence tolerance including WRKY, F-box, zinc fingers, glycosyltransferase, protein kinase, cytochrome P450, PP2C, hypoxia-responsive family, and DUF domain. By haplotype analysis, the 21 DECGs demonstrated distinct genetic differentiation and substantial genetic distance mainly between indica and japonica subspecies. Further, the MQTL7.1 was successfully validated using flanked marker S2329 on a set of genotypes with phenotypic variation. CONCLUSION: This study provides a new perspective on understanding the genetic basis of submergence tolerance in rice. The identified MQTL and novel candidate genes lay the foundation for marker-assisted breeding/engineering of flooding-tolerant cultivars conducive to direct seeding.

Genome-wide association study and development of molecular markers for yield and quality traits in peanut (Arachis hypogaea L.).

BACKGROUND: This study aims to decipher the genetic basis governing yield components and quality attributes of peanuts, a critical aspect for advancing molecular breeding techniques. Integrating genotype re-sequencing and phenotypic evaluations of seven yield components and two grain quality traits across four distinct environments allowed for the execution of a genome-wide association study (GWAS). RESULTS: The nine phenotypic traits were all continuous and followed a normal distribution. The broad heritability ranged from 88.09 to 98.08%, and the genotype-environment interaction effects were all significant. There was a highly significant negative correlation between protein content (PC) and oil content (OC). The 10× genome re-sequencing of 199 peanut accessions yielded a total of 631,988 high-quality single nucleotide polymorphisms (SNPs), with 374 significant SNP loci identified in association with the nine traits of interest. Notably, 66 of these pertinent SNPs were detected in multiple environments, and 48 of them were linked to multiple traits of interest. Five loci situated on chromosome 16 (Chr16) exhibited pleiotropic effects on yield traits, accounting for 17.64-32.61% of the observed phenotypic variation. Two loci on Chr08 were found to be strongly associated with protein and oil contents, accounting for 12.86% and 14.06% of their respective phenotypic variations, respectively. Linkage disequilibrium (LD) block analysis of these seven loci unraveled five nonsynonymous variants, leading to the identification of one yield-related candidate gene and two quality-related candidate genes. The correlation between phenotypic variation and SNP loci in these candidate genes was validated by Kompetitive allele-specific PCR (KASP) marker analysis. CONCLUSIONS: Overall, molecular markers were developed for genetic loci associated with yield and quality traits through a GWAS investigation of 199 peanut accessions across four distinct environments. These molecular tools can aid in the development of desirable peanut germplasm with an equilibrium of yield and quality through marker-assisted breeding.

A membrane associated tandem kinase from wild emmer wheat confers broad-spectrum resistance to powdery mildew.

Crop wild relatives offer natural variations of disease resistance for crop improvement. Here, we report the isolation of broad-spectrum powdery mildew resistance gene Pm36, originated from wild emmer wheat, that encodes a tandem kinase with a transmembrane domain (WTK7-TM) through the combination of map-based cloning, PacBio SMRT long-read genome sequencing, mutagenesis, and transformation. Mutagenesis assay reveals that the two kinase domains and the transmembrane domain of WTK7-TM are critical for the powdery mildew resistance function. Consistently, in vitro phosphorylation assay shows that two kinase domains are indispensable for the kinase activity of WTK7-TM. Haplotype analysis uncovers that Pm36 is an orphan gene only present in a few wild emmer wheat, indicating its single ancient origin and potential contribution to the current wheat gene pool. Overall, our findings not only provide a powdery mildew resistance gene with great potential in wheat breeding but also sheds light into the mechanism underlying broad-spectrum resistance.

Cas9-directed long-read sequencing to resolve optical genome mapping findings in leukemia diagnostics.

Leukemias are genetically heterogeneous and diagnostics therefore includes various standard-of-care (SOC) techniques, including karyotyping, SNP-array and FISH. Optical genome mapping (OGM) may replace these as it detects different types of structural aberrations simultaneously and additionally detects much smaller aberrations (500 bp vs 5-10 Mb with karyotyping). However, its resolution may still be too low to define clinical relevance of aberrations when they are located between two OGM labels or when labels are not distinct enough. Here, we test the potential of Cas9-directed long-read sequencing (LRS) as an additional technique to resolve such potentially relevant new findings. From an internal Bionano implementation study we selected ten OGM calls that could not be validated with SOC methods. Per variant we designed crRNAs for Cas9 enrichment, prepared libraries and sequenced them on a MinION/GridION device. We could confirm all aberrations and, importantly, the actual breakpoints of the OGM calls were located between 0.2 and 5.5 kb of the OGM-estimated breakpoints, confirming the high reliability of OGM. Furthermore, we show examples of redefinition of aberrations between labels that enable judgment of clinical relevance. Our results suggest that Cas9-directed LRS can be a relevant and flexible secondary technique in diagnostic workflows including OGM.

Genetic insights into superior grain number traits: a QTL analysis of wheat-Agropyron cristatum derivative pubing3228.

BACKGROUND: Agropyron cristatum (L.) is a valuable genetic resource for expanding the genetic diversity of common wheat. Pubing3228, a novel wheat-A. cristatum hybrid germplasm, exhibits several desirable agricultural traits, including high grain number per spike (GNS). Understanding the genetic architecture of GNS in Pubing3228 is crucial for enhancing wheat yield. This study aims to analyze the specific genetic regions and alleles associated with high GNS in Pubing3228. METHODS: The study employed a recombination inbred line (RIL) population derived from a cross between Pubing3228 and Jing4839 to investigate the genetic regions and alleles linked to high GNS. Quantitative Trait Loci (QTL) analysis and candidate gene investigation were utilized to explore these traits. RESULTS: A total of 40 QTLs associated with GNS were identified across 16 chromosomes, accounting for 4.25-17.17% of the total phenotypic variation. Five QTLs (QGns.wa-1D, QGns.wa-5 A, QGns.wa-7Da.1, QGns.wa-7Da.2 and QGns.wa-7Da.3) accounter for over 10% of the phenotypic variation in at least two environments. Furthermore, 94.67% of the GNS QTL with positive effects originated from Pubing3228. Candidate gene analysis of stable QTLs identified 11 candidate genes for GNS, including a senescence-associated protein gene (TraesCS7D01G148000) linked to the most significant SNP (AX-108,748,734) on chromosome 7D, potentially involved in reallocating nutrients from senescing tissues to developing seeds. CONCLUSION: This study provides new insights into the genetic mechanisms underlying high GNS in Pubing3228, offering valuable resources for marker-assisted selection in wheat breeding to enhance yield.

PySmooth: a Python tool for the removal and correction of genotyping errors.

In genetic mapping studies involving many individuals, genome-wide markers such as single nucleotide polymorphisms (SNPs) can be detected using different methods. However, it comes with some errors. Some SNPs associated with diseases can be in regions encoding long noncoding RNAs (lncRNAs). Therefore, identifying the errors in genotype file and correcting them is crucial for accurate genetic mapping studies. We develop a Python tool called PySmooth, that offers an easy-to-use command line interface for the removal and correction of genotyping errors. PySmooth uses the approach of a previous tool called SMOOTH with some modifications. It inputs a genotype file, detects errors and corrects them. PySmooth provides additional features such as imputing missing data, better user-friendly usage, generates summary and visualization files, has flexible parameters, and handles more genotype codes. AVAILABILITY AND IMPLEMENTATION: PySmooth is available at https://github.com/lncRNAAddict/PySmooth .

Many purported pseudogenes in bacterial genomes are bona fide genes.

BACKGROUND: Microbial genomes are largely comprised of protein coding sequences, yet some genomes contain many pseudogenes caused by frameshifts or internal stop codons. These pseudogenes are believed to result from gene degradation during evolution but could also be technical artifacts of genome sequencing or assembly. RESULTS: Using a combination of observational and experimental data, we show that many putative pseudogenes are attributable to errors that are incorporated into genomes during assembly. Within 126,564 publicly available genomes, we observed that nearly identical genomes often substantially differed in pseudogene counts. Causal inference implicated assembler, sequencing platform, and coverage as likely causative factors. Reassembly of genomes from raw reads confirmed that each variable affects the number of putative pseudogenes in an assembly. Furthermore, simulated sequencing reads corroborated our observations that the quality and quantity of raw data can significantly impact the number of pseudogenes in an assembler dependent fashion. The number of unexpected pseudogenes due to internal stops was highly correlated (R2 = 0.96) with average nucleotide identity to the ground truth genome, implying relative pseudogene counts can be used as a proxy for overall assembly correctness. Applying our method to assemblies in RefSeq resulted in rejection of 3.6% of assemblies due to significantly elevated pseudogene counts. Reassembly from real reads obtained from high coverage genomes showed considerable variability in spurious pseudogenes beyond that observed with simulated reads, reinforcing the finding that high coverage is necessary to mitigate assembly errors. CONCLUSIONS: Collectively, these results demonstrate that many pseudogenes in microbial genome assemblies are actually genes. Our results suggest that high read coverage is required for correct assembly and indicate an inflated number of pseudogenes due to internal stops is indicative of poor overall assembly quality.

Fine mapping and breeding application of two brown planthopper resistance genes derived from landrace rice.

The Brown planthopper (Nilaparvata lugens Stål; BPH) is known to cause significant damage to rice crops in Asia, and the use of host-resistant varieties is an effective and environmentally friendly approach for controlling BPH. However, genes limited resistance genes that are used in insect-resistant rice breeding programs, and landrace rice varieties are materials resources that carry rich and versatile genes for BPH resistance. Two landrace indica rice accessions, CL45 and CL48, are highly resistant to BPH and show obvious antibiosis against BPH. A novel resistance locus linked to markers 12M16.983 and 12M19.042 was identified, mapped to chromosome 12 in CL45, and designated Bph46. It was finely mapped to an interval of 480 kb and Gene 3 may be the resistance gene. Another resistance locus linked to markers RM26567 and 11MA104 was identified and mapped to chromosome 11 in CL48 and designated qBph11.3 according to the nominating rule. It was finely mapped to an interval of 145 kb, and LOC_Os11g29090 and LOC_Os11g29110 may be the resistance genes. Moreover, two markers, 12M16.983 and 11MA104, were developed for CL45 and CL48, respectively, using marker-assisted selection (MAS) and were confirmed by backcrossing individuals and phenotypic detection. Interestingly, we found that the black glume color is closely linked to the BPH resistance gene in CL48 and can effectively assist in the identification of positive individuals for breeding. Finally, several near-isogenic lines with a 9311 or KW genetic background, as well as pyramid lines with two resistance parents, were developed using MAS and exhibited significantly high resistance against BPHs.

High-quality chromosome-level genome assembly of Nicotiana benthamiana.

Nicotiana benthamiana is a fundamental model organism in plant research. Recent advancements in genomic sequencing have revealed significant intraspecific genetic variations. This study addresses the pressing need for a precise genome sequence specific to its geographic origin by presenting a comprehensive genome assembly of the N. benthamiana LAB strain from the Republic of Korea (NbKLAB). We compare this assembly with the widely used NbLAB360 strain, shedding light on essential genomic differences between them. The outcome is a high-quality, chromosome-level genome assembly comprising 19 chromosomes, spanning 2,762 Mb, with an N50 of 142.6 Mb. Comparative analyses revealed notable variations, including 46,215 protein-coding genes, with an impressive 99.5% BUSCO completeness score. Furthermore, the NbKLAB assembly substantially improved the QV from 33% for NbLAB360 to 49%. This refined chromosomal genome assembly for N. benthamiana, in conjunction with comparative insights, provides a valuable resource for genomics research and molecular biology. This accomplishment forms a strong foundation for in-depth exploration into the intricacies of plant genetics and genomics, improved precision, and a comparative framework.

High-quality faba bean reference transcripts generated using PacBio and Illumina RNA-seq data.

The genome of faba bean was first published in 2023. To promote future molecular breeding studies, we improved the quality of the faba genome based on high-density genetic maps and the Illumina and Pacbio RNA-seq datasets. Two high-density genetic maps were used to conduct the scaffold ordering and orientation of faba bean, culminating in an increased length (i.e., 14.28 Mbp) of chromosomes and a decrease in the number of scaffolds by 45. In gene model mining and optimisation, the PacBio and Illumina RNA-seq datasets from 37 samples allowed for the identification and correction 121,606 transcripts, and the data facilitated a prediction of 15,640 alternative splicing events, 2,148 lncRNAs, and 1,752 fusion transcripts, thus allowing for a clearer understanding of the gene structures underlying the faba genome. Moreover, a total of 38,850 new genes including 56,188 transcripts were identified compared with the reference genome. Finally, the genetic data of the reference genome was integrated and a comprehensive and complete faba bean transcriptome sequence of 103,267 transcripts derived from 54,753 uni-genes was formed.

Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles.

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

High-throughput diagnostic markers for foliar fungal disease resistance and high oleic acid content in groundnut.

BACKGROUND: Foliar diseases namely late leaf spot (LLS) and leaf rust (LR) reduce yield and deteriorate fodder quality in groundnut. Also the high oleic acid content has emerged as one of the most important traits for industries and consumers due to its increased shelf life and health benefits. RESULTS: Genetic mapping combined with pooled sequencing approaches identified candidate resistance genes (LLSR1 and LLSR2 for LLS and LR1 for LR) for both foliar fungal diseases. The LLS-A02 locus housed LLSR1 gene for LLS resistance, while, LLS-A03 housed LLSR2 and LR1 genes for LLS and LR resistance, respectively. A total of 49 KASPs markers were developed from the genomic regions of important disease resistance genes, such as NBS-LRR, purple acid phosphatase, pentatricopeptide repeat-containing protein, and serine/threonine-protein phosphatase. Among the 49 KASP markers, 41 KASPs were validated successfully on a validation panel of contrasting germplasm and breeding lines. Of the 41 validated KASPs, 39 KASPs were designed for rust and LLS resistance, while two KASPs were developed using fatty acid desaturase (FAD) genes to control high oleic acid levels. These validated KASP markers have been extensively used by various groundnut breeding programs across the world which led to development of thousands of advanced breeding lines and few of them also released for commercial cultivation. CONCLUSION: In this study, high-throughput and cost-effective KASP assays were developed, validated and successfully deployed to improve the resistance against foliar fungal diseases and oleic acid in groundnut. So far deployment of allele-specific and KASP diagnostic markers facilitated development and release of two rust- and LLS-resistant varieties and five high-oleic acid groundnut varieties in India. These validated markers provide opportunities for routine deployment in groundnut breeding programs.

A cluster of putative resistance genes is associated with a dominant resistance to sunflower broomrape.

KEY MESSAGE: The HaOr5 resistance gene is located in a large genomic insertion containing putative resistance genes and provides resistance to O. cumana, preventing successful connection to the sunflower root vascular system. Orobanche cumana (sunflower broomrape) is a parasitic plant that is part of the Orobanchaceae family and specifically infests sunflower crops. This weed is an obligate parasitic plant that does not carry out photosynthetic activity or develop roots and is fully dependent on its host for its development. It produces thousands of dust-like seeds per plant. It possesses a high spreading ability and has been shown to quickly overcome resistance genes successively introduced by selection in cultivated sunflower varieties. The first part of its life cycle occurs underground. The connection to the sunflower vascular system is essential for parasitic plant survival and development. The HaOr5 gene provides resistance to sunflower broomrape race E by preventing the connection of O. cumana to the root vascular system. We mapped a single position of the HaOr5 gene by quantitative trait locus mapping using two segregating populations. The same location of the HaOr5 gene was identified by genome-wide association. Using a large population of thousands of F2 plants, we restricted the location of the HaOr5 gene to a genomic region of 193 kb. By sequencing the whole genome of the resistant line harboring the major resistance gene HaOr5, we identified a large insertion of a complex genomic region containing a cluster of putative resistance genes.

A fast comparative genome browser for diverse bacteria and archaea.

Genome sequencing has revealed an incredible diversity of bacteria and archaea, but there are no fast and convenient tools for browsing across these genomes. It is cumbersome to view the prevalence of homologs for a protein of interest, or the gene neighborhoods of those homologs, across the diversity of the prokaryotes. We developed a web-based tool, fast.genomics, that uses two strategies to support fast browsing across the diversity of prokaryotes. First, the database of genomes is split up. The main database contains one representative from each of the 6,377 genera that have a high-quality genome, and additional databases for each taxonomic order contain up to 10 representatives of each species. Second, homologs of proteins of interest are identified quickly by using accelerated searches, usually in a few seconds. Once homologs are identified, fast.genomics can quickly show their prevalence across taxa, view their neighboring genes, or compare the prevalence of two different proteins. Fast.genomics is available at https://fast.genomics.lbl.gov.

Genetic linkage mapping and quantitative trait locus (QTL) analysis of sweet basil (Ocimum basilicum L.) to identify genomic regions associated with cold tolerance and major volatiles.

Chilling sensitivity is one of the greatest challenges affecting the marketability and profitability of sweet basil (Ocimum basilicum L.) in the US and worldwide. Currently, there are no sweet basils commercially available with significant chilling tolerance and traditional aroma profiles. This study was conducted to identify quantitative trait loci (QTLs) responsible for chilling tolerance and aroma compounds in a biparental mapping population, including the Rutgers advanced breeding line that served as a chilling tolerant parent, 'CB15', the chilling sensitive parent, 'Rutgers Obsession DMR' and 200 F2 individuals. Chilling tolerance was assessed by percent necrosis using machine learning and aroma profiling was evaluated using gas chromatography (GC) mass spectrometry (MS). Single nucleotide polymorphism (SNP) markers were generated from genomic sequences derived from double digestion restriction-site associated DNA sequencing (ddRADseq) and converted to genotype data using a reference genome alignment. A genetic linkage map was constructed and five statistically significant QTLs were identified in response to chilling temperatures with possible interactions between QTLs. The QTL on LG24 (qCH24) demonstrated the largest effect for chilling response and was significant in all three replicates. No QTLs were identified for linalool, as the population did not segregate sufficiently to detect this trait. Two significant QTLs were identified for estragole (also known as methyl chavicol) with only qEST1 on LG1 being significant in the multiple-QTL model (MQM). QEUC26 was identified as a significant QTL for eucalyptol (also known as 1,8-cineole) on LG26. These QTLs may represent key mechanisms for chilling tolerance and aroma in basil, providing critical knowledge for future investigation of these phenotypic traits and molecular breeding.

Genomic basis determining root system architecture in maize.

KEY MESSAGE: A total of 389 and 344 QTLs were identified by GWAS and QTL mapping explaining accumulatively 32.2-65.0% and 23.7-63.4% of phenotypic variation for 14 shoot-borne root traits using more than 1300 individuals across multiple field trails. Efficient nutrient and water acquisition from soils depends on the root system architecture (RSA). However, the genetic determinants underlying RSA in maize remain largely unexplored. In this study, we conducted a comprehensive genetic analysis for 14 shoot-borne root traits using 513 inbred lines and 800 individuals from four recombinant inbred line (RIL) populations at the mature stage across multiple field trails. Our analysis revealed substantial phenotypic variation for these 14 root traits, with a total of 389 and 344 QTLs identified through genome-wide association analysis (GWAS) and linkage analysis, respectively. These QTLs collectively explained 32.2-65.0% and 23.7-63.4% of the trait variation within each population. Several a priori candidate genes involved in auxin and cytokinin signaling pathways, such as IAA26, ARF2, LBD37 and CKX3, were found to co-localize with these loci. In addition, a total of 69 transcription factors (TFs) from 27 TF families (MYB, NAC, bZIP, bHLH and WRKY) were found for shoot-borne root traits. A total of 19 genes including PIN3, LBD15, IAA32, IAA38 and ARR12 and 19 GWAS signals were overlapped with selective sweeps. Further, significant additive effects were found for root traits, and pyramiding the favorable alleles could enhance maize root development. These findings could contribute to understand the genetic basis of root development and evolution, and provided an important genetic resource for the genetic improvement of root traits in maize.